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Abstract 

Recently, several theoretical and experimental studies have been undertaken to 
probe the effect of stochasticity on gene expression (GE). In experiments, the GE re- 
sponse to an inducing signal in a cell, measured by the amount of mRN As /proteins 
synthesized, is found to be either graded or binary The latter type of response 
gives rise to a bimodal distribution in protein levels in an ensemble of cells. One 
possible origin of binary response is cellular bistability achieved through positive 
feedback or autoregulation. In this paper, we study a simple, stochastic model of 
GE and show that the origin of binary response lies exclusively in stochasticity. 
The transitions between the active and inactive states of the gene are random in 
nature. Graded and binary responses occur in the model depending on the rel- 
ative stability of the activated and deactivated gene states with respect to that 
of mRN As /proteins. The theoretical results on binary response provide a good 
description of the "all-or-none" phenomenon observed in an eukaryotic system. 

Keywords: gene expression, graded and binary responses, stochastic binary 
response, "all-or-none" phenomenon, probability density, activation 

1. Introduction 

Gene expression (GE) is the central activity in a living cell. The two ma- 
jor steps in GE, transcription and translation, involve several biochemical 
reactions. The time evolution of this system of reactions or events is not a 
continuous process as molecular population levels in a reacting system change 
only by discrete amounts. Furthermore, the time evolution is not determin- 
istic as the biochemical events underlying GE are probabilistic in nature, 
i.e., the timing of the biochemical events cannot be predicted with certainty. 
For example, the binding/unbinding of RNA polymerase (RNAP) at the 
promoter region of DNA and that of regulatory molecules at the operator 
regions are probabilistic processes. The discrete, probabilistic nature of the 



biochemical events may be ignored in the limit of large numbers of partici- 
pating biochemical molecules. In this case, the biochemical reactions/events 
occur at much higher frequencies and fluctuations around the mean levels of 
biomolecules participating in GE are small. Thus, the time evolution of the 
system of reactions may be treated to be continuous and deterministic as in 
the traditional differential rate equation approach. In a living cell, the num- 
ber of biomolecules involved in GE is often small so that a stochastic rather 
than deterministic description provides the more correct picture. In recent 
years, there is a growing realization that stochasticity plays an important 
role in determining the outcome of biochemical processes in the cell [TJ [2]. 
Stochastic effects in GE explain the pronounced cell-cell variation observed in 
isogenic populations. A cell may have the option of proceeding along one of 
two possible developmental pathways. The pathway selection is probabilistic 
and the cell fate depends on the particular choice of pathway. Thus, even a 
clonal population of cells can give rise to two distinct subpopulations in the 
course of time. The randomization of pathway choice leads to diversity and 
increases the likelihood of survival of organisms in widely different environ- 
ments. A well-known example of the two way-choice is that of lysis-lysogeny 
in E. Coli 0. 

The effect of stochasticity (randomness/noise) is prominent at the level 
of an individual cell and can be masked due to ensemble averaging in a 
population of cells. Single cell experiments provide evidence that GE in a 
cell occurs in abrupt stochastic bursts (HEJEI- In more recent experiments, 
a quantitative measure of noise associated with GE has been obtained in 
both prokaryotic as well as eukaryotic cells HI EH A large number of 
theoretical studies address the origin and consequences of stochasticity in 

GE |Sl HOI HH H21 C31 IHl ESI EE! EH CHI - Thus ; the notion of stochasticity in 
GE is well established both theoretically and experimentally. 

Regulation of GE in a cell is achieved in a manifold of ways which in- 
crease in complexity from prokaryotic to eukaryotic cells. In the prokaryotic 
systems, regulation is achieved by the binding of regulatory molecules (re- 
pressors or activators) to the operator regions of DNA. In eukaryotes, the 
activator molecules are known as transcription factors (TFs). Intra- and 
extra- cellular inducing signals activate the TFs which then bind to appro- 
priate enhancer sequences on the DNA. The GE response to an inducing 
signal in an individual cell may be graded or binary. Response is quantified 
by the amount of mRN As / proteins synthesized. In graded response, the 
output varies continuously as the amount of input stimulus is varied till the 
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steady state is reached. In binary response, alternatively termed the "all- 
or-none" phenomenon, the output has a binary character, i.e., GE occurs at 
either a low or a high level and expression at intermediate levels is rare. This 
gives rise to a bimodal distribution in protein levels in an ensemble of cells. 
Several experiments on both prokaryotic and eukaryotic cells establish the 
binary character of GE [BJ 03 EH EDI EH 1221 EH • Binary response may be as- 
cribed to bistability which implies existence of two stable steady states with 
low and high protein levels of GE. One way of achieving bistability is through 
positive feedback or autocatalysis in which the protein product of GE pro- 
motes further GE either directly or via intermediates. The lac operon in E. 
Coli is an example of a model system in which autocatalytic induction gives 
rise to the "all-or-none" phenomenon in GE [EH 123 EH EH1 EH • Beckskei et al. 
[23j have demonstrated that positive feedback can generate binary response 
in a synthetic eukaryotic gene circuit. In eukaryotic transcription, enhancers 
activate the usually weak eukaryotic promoters. There is now strong experi- 
mental evidence that in some systems enhancers do not affect transcription 
rate but rather increase the probability of a gene synthesizing proteins at a 
high level (23 EHl EH EHj- In a population of cells, enhancers increase the 
number of cells expressing at a high level but not the level of expression per 
cell. 

Binary response in GE can have a purely stochastic origin. Kepler and 
Elston (TO] provide examples of stochastic binary response (SBR), i.e., binary 
response induced by noise. A simple model of SBR shows a binary distribu- 
tion of mRNA levels in an ensemble of cells [HI]. A recent model of eukary- 
otic GE suggests that fluctuations in the binding of TFs to DNA can explain 
graded and binary responses observed in inducible GE |32|. Fast chemical 
kinetics is responsible for a graded response whereas slow kinetics leads to a 
binary output. The "all-or-none" phenomenon observed in some eukaryotic 
systems does not involve positive feedback processes explicitly jUIHlE]. On 
the other hand, protein synthesis in these systems occurs in stochastic bursts. 
Since the effect of stochasticity is prominent in these systems, it is reasonable 
to conjecture that the "all-or-none" phenomenon (binary response) observed 
in these systems is a manifestation of stochasticity. In this paper, we con- 
sider a simple, stochastic model of GE studied earlier [T3 EH] • We show 
that graded and binary responses occur naturally in the model depending on 
the relative stability of activated and deactivated gene states with respect to 
that of mRN As /proteins. Binary response, obtained in the model, arises 
solely due to stochasticity and not due to positive feedback processes. We 
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further show that our model gives a good description of the "all-or-none" 
phenomenon observed in an eukaryotic system P|. 

2. Stochastic model of GE 

In the minimal model of GE, a gene can be in two possible states: inactive 
(G) and active (G*). Random transitions occur between the states G and 
G* according to the first order kinetics 



k a jp lip 

G ^ G* — > p — > $ (1) 
k d 

where k a and kd are the activation and deactivation rate constants. In the 
active state G*, transcription is initiated followed by translation and protein 
synthesis. The separate processes are combined into a single step of protein 
(p) synthesis with the rate constant j p . The protein degrades with the rate 
constant k p and the degradation product is represented as $. If cell division 
is taken into account, the protein decay rate has two components, one the 
degradation rate and the other the dilution rate of proteins due to cell growth 
and division. In this case, k p denotes the rate constant for protein decay. 

In inducible GE systems, the activation of a gene is brought about by an 
activator S, say, TFs. The reaction scheme in the presence of S is given by 

k\ k a jp kp 

G + S ^ GS ^ G* — > p — > $ (2) 
k 2 k d 

where GS represents the bound complex of G and S. The reaction scheme 
in equation (|2j) can be generalized by including direct transitions between G 
and G*. The rate constants for transitions from G to G* and G* to G are k on 
and k a ff respectively. For eukaryotic systems, the rate constant k on has a 
very low value as activating TFs, S, are required in most cases for transistion 
to the active state G*. The reaction scheme is given by 

k\ k a k on 

G + S ^ GS ^ G* , G ^ G* 
k 2 k d k ff 
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G* — > p — > $ 1 ; 

For eukaryotic systems, the initiation of transcription by RNA polymerase 
II generally requires a prior assembly of TFs on the enhancer regions of the 
target gene. This state of the gene is represented by G-S in the reaction 
scheme 03 The activating TFs, S, facilitate the formation of the transcription 
initiation complex which is bound to the promoter region of DNA and con- 
sists of general TFs, other factors and RNA polymerase II. The gene is now 
in the active state G* and RNAP starts transcription after disengaging itself 
from the initiation complex through the key step of phosphorylation. The 
general TFs are then released allowing for the initiation of a new round of 
transcription with another RNAP molecule. In the simple reaction scheme 
El this is respresented by a return after transcription initiation to the inter- 
mediate complex G-S and subsequent return to the active state G*. 
If n G be the total concentration of genes then 

n G = [G] + [G-S] + [G*] (4) 

where [G], [G-S], and [G*] denote the concentrations of genes in the states 
G, G-S, and G* respectively. Using the method of King and Altman [34|, 
the fractions of genes in the inactive, intermediate and active states are given 
by 

[G] _ fa k off + fa faff + fa fa 
n G A 
[G-S] _k d k l S + k on fa + k of f fa S 



n G A 

[G*] k a k\ S + k a k on + fa k 

n G A 



(5) 



respectively, where 



A — k a fa S + fa k on + fa k on + k/j, fa S-\- k on k d + k Q ff fa S + fafaff + fa faf f + fafa 

(6) 

From equation (jSJ), one can further write 
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n G k a 



(7) 



where 



Also, 



(r + D + kon , , k a /k' + k of f (1 + f ) 

«. A T & T kj \ ~ k ~ k s > 



ky, &2 ' ^2 , , 

k s = — , re = - — ana k = (9) 

k\ k on ^off 



In the reaction scheme HI the steady state concentration of genes in the active 
state is given by 

Expressions Q and (|T0l) are identical in form with k a and kd replaced by k a , 
k d . The equivalence relations in equation (jBJ) enable one to map the reaction 
scheme El onto the simpler schemed] while calculating various quantities. Use 
of the simpler reaction scheme leads to greater mathematical tractability. 
The half-lives of the active and inactive states of the gene in the reaction 
scheme El are given by T a = log2/k a and T d = log2/k d respectively. Since k a 
and k d are given by equation (jBJ), the half-lives are dependent on k a , kd as 
well as S, the concentration of TFs. 

We now consider a simple stochastic model corresponding to reaction 
scheme [IJ The results we derive hold true for the more complicated reaction 
scheme El but with k a , k d replaced by k a , k d (equation ((HJ)). At this point, one 
can ask about the validity of the equivalence relations (equations (7) and (8)) 
in the stochastic case. Use of the relations is justified only if the fluctuations 
in the concentration S of the activator molecules are ignored. Exact valid- 
ity can be established by deriving expressions for variance from the Master 
Equations (treating S to be constant) corresponding to the reaction schemes 
1 and 3. This has been done for the simpler case k on = 0, k Q ff = (in the 
general case, these rate constants are much smaller than the activation and 
deactivation rate constants k a and k d ). The expressions for variance in the 
reaction schemes 1 and 3, are found to be identical with k a , k d in scheme 1 
replaced by k a , k d in scheme 3. In the model, the only stochasticity arises 
from random transitions of a gene between the inactive and active states as 
in the minimal model of Cook et al. ^H]- Protein synthesis from the active 
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gene and protein degradation occur in a deterministic manner. In each state 
of the gene, the concentration of proteins evolves deterministically according 
to the equation 

— = -y^z - k p x = f(x,z) (11) 

-^-max 

where z = 1(0) when the gene is in the active (inactive) state and x = x , 
X and X max being the protein concentration at time t and the maximum 
protein concentration respectively. The variable x thus denotes protein con- 
centration normalized by the maximum possible concentration. The latter 
quantity is equal to the protein concentration in the steady state if the gene is 
always in the active state, i.e, deactivation processes are disallowed. We note 
that X max = |r. Let pj(x,t) (j = 0, 1) be the probability density function 
when z = j. The total probability density function is 

p(x 1 t)=p (x,t)+pi(x,t) (12) 

The rate of change of probability density is given by 

= -^\fM)Pii.^)] + EtNWM) - W jkPj {x,t)] (13) 

X kjtj 

where Wkj is the transition rate from the state k to the state j and Wjk is 
the same for the reverse transition. The first term in equation (fT3j) is the so 
called "transport" term representing the net flow of the probability density. 
The second term represents the gain/loss in the probability density due to 
random transitions between the state j and other accessible states. In the 
present case, equation (fT3|) gives rise to the following two equations: 

dP °)^ - = -J-(-^PoOM)) + k d p!(x,t) - k a Po(x,t) (14) 



dPl< f' - = ~"J-{(y^ k p x)pi(x,t)} + k aPo (x,t) - kdPi(x,t) (15) 

OX OX -A- max 

The steady state distribution ( ap °^'^ = 0, dpi ^ ,e> = 0) is given by 

p(x) = N x { ^~ 1] (I - x) { ^~ l) (16) 
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(a) (b) (c) 

Figure 1: Plot of p(x) versus x (Case I, r\, r 2 > 1): 1A (r 2 >> r-y), IB 
(r 2 = ri) and 1C (7*1 >> r 2 ) respectively. 



where N, the normalization constant, is given by the inverse of a beta func- 
tion m 

N = -r^r, (i7) 



B( 



3. Graded and binary responses 

The graded and binary responses to varying concentrations of S can be un- 
derstood by considering two limiting cases of the steady state distribution of 
protein levels p(x) (equation ([IB]) ): |s > 1 (Case I) and |* < 1 (Case 
II). Figure 1 shows plots of p(x) versus x corresponding to Case I. The mean 
value of x is given by 

< x > = / x p(x) dx (18) 



k p k a + kd 

Define r\ = I s - and r 2 = -A. Figures 1A, IB and 1C correspond to 

r 2 >> r\ > 1, ri = r 2 > 1 and ri >> r 2 > 1 respectively. In the presence of 

an inducing stimulus, k a and kd are replaced by k a and k d (equation (JHJ)). 

As ri = j5* increases, the mean protein level increases from a lower to a 
higher value. The increase in r\ can be brought about by increasing the 
concentration of S. Thus the mean protein level is a continuous function of 
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[S], i.e., a graded response is obtained. Saturation level is reached when jj- 
in equation (jBJ) is >> 1 so that k' a = k a and k' d = k d + k a ff . 

Figure 2 shows plots of p(x) versus x corresponding to Case II, i.e., t\ < 1, 
r 2 < 1. In this case p(x) is peaked at a low (zero) value of a; (r\ « r 2 , figure 
2A), a high value of x (r\ » r 2 , figure 2D) or at both low and high values 
of x (figures 2B and 2C). Thus, in a cell GE predominantly occurs at low 
and/or high levels and protein production at intermediate levels is negligible. 

Again, in the presence of an inducing stimulus, S, r± = and r 2 = -jf can 
be changed by changing the concentration [S]. The response in this case is 
not graded as the mean protein level is not a continuous function of [S] but 
has only low and high values. Figures 2B and 2C correspond to SBR and 
bifurcation from a unimodal probability distribution function (figure 2A) to 
a bimodal one is brought about by varying r\ and r 2 . SBR gives rise to the 
"all-or-none" phenomenon in GE. In experiments on a population of cells, a 
fraction of cells is found to be in the state with low (zero) level of GE and 
another fraction is in the state with high level of GE. The fraction of cells in 
which protein synthesis occurs at intermediate levels is small. In the cases 
when ti > 1, r% < 1 and r\ < 1, r 2 > 1, unimodal responses are obtained. In 
the first case, GE occurs at a high level and in the second case, GE occurs 
at a low level. The response is not graded in the presence of an inducing 
stimulus. 

The graded and binary responses to an inducing stimulus are a natural 
outcome of stochastic gene activation and deactivation processes. If the gene 
is always in the inactive state [z = in equation (jTTl) ). the mean protein 
level corresponds to x = 0. If the gene is in the active state (z = 1) and no 
deactivation processes are allowed, the mean protein level is given by jr and 
x = 1 corresponding to maximum protein synthesis. When stochastic GE 
is considered, i.e., random activation/deactivation processes are taken into 
account, two possibilities arise. If the activation and deactivation rates are 
faster than the protein degradation rate, an average protein level is obtained 
due to the accumulation of proteins over random transitions between the 
values x = and x — 1. In the opposite case, i.e., when the activation and 
deactivation rates are slower than the protein degradation rate, the mean 
protein level is either x = or x = 1 depending on whether the gene is in 
the inactive or the active state. The half-life of each such state is larger than 
that of synthesized proteins so that in each case sufficient time is available 
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for the mean protein level to attain its particular steady state value. Due 
to the relatively larger protein degradation rate, there is no accumulation 
of proteins over the random transitions so that observed protein levels are 
predominantly at x = and x — 1. 

Ko [TT] has considered a stochastic model for gene induction and has 
shown using computer simulation that two types of response are possible. GE 
in the model is switched on and off due to the binding and unbinding of the 
TF-complex at the gene. Stochasticity is introduced into the model because 
of the probabilistic nature of the binding and unbinding events. An unstable 
transcription complex causes a "homogeneous" level of gene induction while 
a stable transcription complex gives rise to a "heterogeneous" level. The 
homogeneous case is analogous to graded response and the binary response 
is an example of heterogeneous response. In the detailed stochastic model 
studied by Pirone and Elston [32j, fluctuations in TF binding are shown to 
explain graded and binary responses to an inducing stimulus. A binary pat- 
tern of GE is obtained when the enhancer-state fluctuations (caused by the 
binding and unbinding of TFs) are slow whereas faster enhancer-state fluctu- 
ation give rise to a graded response. The conclusions are arrived at by using 
a combination of approximate analytical methods and numerical techniques 
like Monte Carlo (MC) simulation based on the Gillespie Algorithm. The 
role of operator fluctuations in transcriptional regulation has been studied 
by Kepler and Elston (TOJ using the Master Equation Approach. In the limit 
of large protein abundance, an equation similar to equation (THfl) is obtained. 
Again, the interpretation is the same. In each state of the operator, the pro- 
tein concentration evolves deterministically but there are random transitions 
between the two states of the operator corresponding to the occupation and 
unoccupation of the operator region by an activator. The analysis is not, 
however, extended further. 
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Figure 2: Plot of p(x) versus x (Case II, ri, r 2 < 1): 2A (r 2 >> ri), 2B 
(r 2 > ri), 2C (r 2 = ri), and 2D (ri >> r 2 ) respectively. 

Discovery of biologically active molecules, say, drugs involves testing the 
effect of such molecules on appropriate targets. Membrane receptors are the 
largest class of drug targets. Drugs interacting with receptors are broadly 
of two types: antagonists and agonists. Antagonists block receptor activity 
whereas agonists binding to the receptors enhance cellular activity. The bind- 
ing triggers a series of biochemical events which lead to a change in cellular 
activity. The change can be linked to the expression of a reporter gene so 
that detection and quantification of the response to agonist-induced receptor 
activation are possible. Figure 3 shows a cartoon of how the reporter gene 
conveys information regarding receptor activation pj. A cascade of intra- 
cellular processes are initiated by the binding of an agonist to the receptor. 
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Figure 3: Schematic diagram showing the action of carbachol. 

This is accompanied by a change in the concentrations of messengers in the 
cell. As a result, TFs are activated which then translocate to the nucleus. 
The TFs bind to the target gene and initiate expression of both this gene 
as well the reporter gene. The mRNA, generated from the reporter gene, 
is translated into the reporter enzyme. The reporter enzyme catalyzes the 
conversion of substrates into detectable products. 

The scenario depicted in figure 3 provides the basis for high-throughput 
screening of pharmaceutical candidate drugs in living mammalian cells 
The reporter gene used is (3 — lactamase the protein product of which hy- 
drolyzes a substrate giving rise to a large shift in fluorescence emission wave- 
length. Cells in which the reporter gene is not expressed or expressed at a 
very low level appear fluorescent green whereas reporter-positive cells with a 
high level of GE appear fluorescent blue. The activation of cellular processes 
is brought about by the binding of the agonist carbachol to the muscarinic 
receptor. In the experiment, the percentages of blue, blue-green (intermedi- 
ate level of GE) and green cells are measured by flow cytometry with varying 
carbachol dose and also as a function of time after stimulation by carba- 
chol. The major finding is that as the carbachol dose increases from a low 
to high value, the fraction of green cells (low level of GE) decreases and that 
of blue cells (high level of GE) increases. The percentage of blue-green cells 
remains fairly low throughout. This is a manifestation of the "all-or-none" 
phenomenon, i.e., binary response in GE. 
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We now show that the simple stochastic model studied by us provides a 
good quantitative fit to the experimental data of Zlokarnik et al. Since the 
probability density function p(x) of protein levels is known (equation JUJ), 
one can calculate experimentally measurable quantities like the dose-response 
function. Figure 4 shows the experimental data points corresponding to 
fractions of blue + blue-green (depicted by solid triangles) and blue (depicted 
by solid squares) cells versus log (C), where C is the carbachol concentration. 
The fraction of blue-green cells is given by the difference in data points 
belonging to the two curves. The remaining cell fractions describe green 
cells. The concentration of activated TFs (S in our model) may be taken to 
be proportional to the concentration C of carbachol and in our theoretical 
dose-response curves (solid lines in figure 4), [S] replaces C. From equation 
(fT6|) . the steady state probability of finding a cell with x (mean protein level 
divided by maximum protein concentration), greater than a threshold value 

Xthr IS 

[* thT (l-x) ( % _1) dz , s 

p(x > x thr ) = 1 - ^ - h 20 

J x k p (1 - x) k p dx 

^ P thr k p ' k p ' k p ' tnr > (21) 

kaB(t p , ft) 

where 2-^1 ( a > b, c; z) is the hypergeometric function [35] . In our model, we 
assume that a cell is in a state with high level of GE if the mean protein 
level in the steady state is greater than a fraction of 0.9 of the maximum 
protein concentration i.e., x > 0.9. By setting x t hr = 0.9 in equation (|2T)1) 
and replacing k a , k d by the effective rate constants k a , k d (equation ©), one 
can calculate p(x > x t hr) for various values of S. The probability p(x > x t hr) 
can also be interpreted as the fraction of cells in a cell population with 
x > Xthr- The theoretical dose-response curve obtained in this manner gives 
a good fit to the experimental data points (solid squares in figure 4) for the 
parameter values (in arbitrary units) k 2 = 1.6 x 10~ 4 , k on = 1.2 x 10~ 6 , 
k off = 1.32 x 10~ 4 , k s = 1.6 x 10- 6 , k p = 1, k a = 0.17 and k d = 0.0465. The 
data points in this case correspond to the fraction of cells in a high level of 
GE (blue cells). A cell is assumed to be in a state with low level of GE if x 
is < x t hr = 0.1 (green cells). A cell is in a state with intermediate level of 
GE when 0.1 < x < 0.9 (blue-green cells). The cell fraction in the last case 
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is given by 



p(0.1 < x < 0.9) 



J 0.1 



-1) 



-i) 



fix 



/ fcq 
X fc P 



-1) 



1 



1) 



c/.r 



(22) 



with fc a , fed replaced by k' a and fc d . The fraction of blue + blue-green cells is 
computed from an expression similar to 22 but with the integration limits 
(0.1, 0.9) in the numerator replaced by (0.1, 1.0). The calculated dose- 
response curve gives a good fit to the experimental data points (solid trian- 
gles in figure 4). The two curves in figure 4 have been obtained for the same 
set of parameter values using Mathematica. The good quantitative agree- 
ment between our theoretical results and experimental data indicates that 
the stochastic model of GE considered by us captures the essential features 
of stochastically induced binary response in GE. 
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Figure 4: Fractions of blue + blue-green (upper curve) and blue (lower curve) 
cells versus S in a semi-logarithm plot. The experimental points are depicted 
by solid triangles signs and solid squares. The parameter values used for the 
fitting are mentioned in the text. 

The reaction scheme 3 doesnot fulfill detailed balance, i.e, equilibrium 
conditions. The steady state concentrations in equation (5) are derived from 
the more general non-equilibrium conditions of stationarity. Ref. [SH] pro- 
vides examples of reaction schemes related to GE which violate detailed 
balance conditions. One of these involves the phenomenon of stochastic fo- 
cussing (SF) in which signal fluctuations sharpen the response in a regulated 
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process. SF appears to be an out-of-equilibrium effect which is absent if the 
reaction scheme is constrained by detailed balance. Modifications of reac- 
tion scheme 3 which preserve detailed balance are possible and it will be 
of interest to determine whether SBR can occur in both equilibrium and 
out-of-equilibrium scenarios. 

4. Conclusion and outlook 

In this paper, we have considered a simple stochastic model of GE and demon- 
strated that stochasticity provides the basis for graded and binary responses 
to inducing signals. The sole ingredients of the minimal model of stochastic 
GE studied in this paper are gene activation, deactivation, protein synthe- 
sis and degradation, each of which involves a number of biochemical events. 
Stochasticity in this model is associated only with the gene activation and de- 
activation processes whereas protein synthesis and degradation are assumed 
to occur in a deterministic manner. A deterministic description of protein 
synthesis is justified when the number of proteins produced is large. This is 
the situation in the experiment by Zlokarnik et al. [H] in which proteins per 
cell are a few thousands in number. For smaller protein numbers, the inclu- 
sion of stochasticity during protein synthesis and degradation is expected to 
blur the GE responses but the major conclusions of the paper still remain 
valid. The processes of transcription and translation in the model are not 
treated separately but lumped together in a single protein synthesis step. 
In an eukaryotic cell, combining transcription and translation into a single 
step may be considered to be a drastic approximation. One can study the 
effect of stochastic gene activation and deactivation on the transcription pro- 
cess itself and focus on mRNA synthesis rather than proteins in reaction 
schemes [TJOl This type of approach highlights the quantal nature of tran- 
scription with bursts of mRNAs being produced in a probabilistic manner in 
agreement with experimental observations [HUSH! . In fact, the value k p = 1 
is more appropriate if k p is interpreted as the mRNA, rather than protein 
decay constant. The mathematical analysis and conclusions are the same as 
before since protein production is linked to mRNA synthesis. Despite the 
limitations of the model, it contains the important features necessary for an 
explanation of the stochastically induced "all-or-none" phenomenon observed 
in some eukaryotic systems. The model results give a good description of 
the experimental data of Zlokarnik et al. || and are expected to be of rele- 
vance in explaining the binary response in GE observed in other eukaryotic 
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systems [HE]. The probabilistic nature of gene activation and deactivation 
processes is crucial to explain how graded and binary responses in GE occur 
in the model. The stochastic origin of binary response is distinctive from the 
binary response brought about by positive feedback processes. Experiments 
designed to probe the stochastic origins of graded and binary responsPlot 
of p(x) versus x (Case I, r 1 , r 2 > 1): 1A (r 2 >> ri), IB (r 2 = r\) and 1C 
(7*1 >> r 2 ) respectively.es, are needed for a clearer understanding of the role 
of stochasticity in such responses. The stochastic model of GE, corresponding 
to reaction scheme 2, has earlier been studied to explore the stochastic ori- 
gins of haploinsufficiency [THl [T3 Qlj. The model studied in the paper is 
a modification of the earlier model. The simplicity of the models allows 
for mathematical analysis and helps in identifying the origins of phenom- 
ena associated with stochastic GE. The knowledge and insight gained from 
the study of simple models like the present one provide necessary inputs to 
develop more detailed and realistic models of GE. 
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